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Abstract 

During the course of our optical deep survey program on La emitters at z ~ 5.7 
in the sky area surrounding the quasar SDSSp J104433.04—012502.2 at z = 5.74, we 
found that a faint galaxy with ms{ AB) ~ 25 is located at l."9 southwest of the quasar. 

Its broad-band color properties from B to z' suggest that the galaxy is located at a 
redshift of z ~ 1.5 — 2.5. This is consistent with no strong emission line in our 
optical spectroscopy. Since the counter image of the quasar cannot be seen in our 
deep optical images, the magnification factor seems not to be very high. Our modest 
estimate is that this quasar is gravitationally magnified by a factor of 2. 

Key words: gravitational lensing — galaxies: high-redshift — galaxies: quasars: 
individual (SDSSp J104433.04-012502.2) 


Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory 
of Japan. 
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1. Introduction 


The Sloan Digital Sky Survey (SDSS: e.g., York et al. 2000) has been finding high- 
redshift quasars at z ~ 6; the most distant known to date is SDSSp J103027.10+052455.0 at 
z = 6.28 (Fan et al. 2000, 2001). One serious problem seems to be that all of these z ~ 6 
SDSS quasars are exceptionally bright (see table 1). As claimed by Wyithe and Loeb (2002; 
hereafter WL02), the presence of such very bright quasars at z ~ 6 raises the following two 
problems: (1) If the central engine of these quasars is an accreting supermassive black hole 
and shines near the Eddington accretion rate, a supermassive black hole with mass exceeding 
~ 3 x 1O 9 M 0 was already formed beyond z — 6. This challenges models for early structure 
formation (Turner 1991; Haiman, Loeb 2001). Also, (2) the absolute luminosities of the z ~ 6 
quasars are systematically higher than expected from the SDSS survey criteria (Fan et al. 2000, 
2001), suggesting that the luminosity function of the z ~ 6 quasars is significantly biased to 
higher luminosities. These two problems urge WL02 to propose that the SDSS z ~ 6 quasars 
may be magnified by intervening gravitational lenses. 

During the course of our optical deep survey program on La emitters at z ~ 5.7 in the 
field surrounding the quasar SDSSp J104433.04—012502.2 at z — 5.74, we found that a faint 
galaxy with m#(AB) ~ 25 is located at l."9 southwest of this quasar. In this paper, we report on 
the possibility that SDSSp J104433.04—012502.2 is gravitationally magnified by a faint galaxy 
with optical photometric and spectroscopic data obtained by our observations. Throughout 
this paper, we adopt a flat universe with D m = 0.3, Da = 0.7, and h = 0.7, where h = Hq /(100 
km s -1 Mpc -1 ). 

2. Optical Deep Imaging 

Deep optical imaging observations were made with the Suprime-Cam (Miyazaki et al. 
1998) on the 8.2 m Subaru telescope (Kaifu 1998) at Mauna Kea Observatories. The Suprime- 
Cam consists of ten 2kx4k CCD chips and provides a very wide field of view: 34' x 27' 
with a Of'2/pixel resolution. Therefore, the combination between the Suprime-Cam and the 
Subaru telescope enables us to carry out wide and deep narrow-band imaging surveys for high- 
z emission-line objects. Using this facility, we made a very deep optical imaging survey for 
faint La emitters in the field surrounding the SDSSp J104433.04—012502.2 at a redshift of 5.74 
(Fan et al. 2000; Djorgovski et al. 2001; Goodrich et al. 2001). In this survey, we used a 
narrow-passband filter, NB816, centered at 8160 A with a passband of AA(FWHM) = 120 
A; the central wavelength corresponds to a redshift of 5.72 for La emission. We also used 
broad-passband filters: B , Rq, Iq, and z'. A summary of the imaging observations is given in 
table 2. All of the observations were made under photometric conditions, and the seeing size 
was between 0''7 and l."3 during the run. 

In this paper, we analyze only two CCD chips in which the quasar SDSSp 
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J104433.04—012502.2 is included (see also Ajiki et al. 2002). The field size is 1U67 by 11.'67. 
The individual CCD data were reduced and combined using IRAF and mosaic-CCD data- 
reduction software developed by Yagi et al. (2002). The photometric and spectrophotometric 
standard stars used in the flux calibration were SA101 for the B, Rc, and Ic data, and GD 
50, GD 108 (Oke 1990), and PG 1034+001 (Massey et al. 1996) for the NB816 data. The z' 
data were calibrated using the magnitude of SDSSp J104433.04—012502.2 (Fan et al. 2000). 
Thumb-nail optical images of the quasar filed are shown in figure 1. A faint galaxy is seen at 
lf'9 southwest from the quasar. This galaxy is seen even in the R-band image, giving a rough 
constraint such that its redshift is z < 2.6, since the wavelength of redshifted L a absorptions 
by intergalactic neutral hydrogen is considered to be shorter than about 4400 A. 

3. Optical Spectroscopy 

In order to estimate a spectroscopic redshift, we carried out optical spectroscopy using 
the Subaru Faint Object Camera And Spectrograph (FOCAS; Kashikawa et al. 2000) on 2002 
March 13 (UT). We used 300 lines mm' 1 grating blazed at 7500 A (300R) together with an 
order cut filter SY47, giving a spectral resolution of R ~ 1000 with a 0."8-wide slit. This setting 
gave a wavelength coverage between 4700 A and 9400 A. The integration time was 1800 s. 
The spectrum is shown in figure 2. This spectrum shows that this galaxy has blue colors, and 
that there is no prominent emission line. If this galaxy is a strong emission-line galaxy, like a 
starburst galaxy, it is inferred either that its redshifted wavelength of La emission line would 
be shorter than 4700 A, or that its redshifted wavelength of [O ll]A3727 emission line would be 
longer than 9400 A. From these constraints, one may estimate that the redshift of this galaxy 
lies in a range between z « 1.5 and z ~ 2.9. 

4. Results and Discussion 

The faint galaxy which we have found (hereafter, the lens galaxy) is located at l."9 
southwest from the quasar. In order to estimate how the quasar is gravitationally magnified, 
we need information about both the redshift and the stellar velocity dispersion of this lens 
galaxy. However, no direct information on these two quantities was obtained in our observations. 
Therefore, in this paper, first, we estimate its photometric redshift based on our optical broad¬ 
band photometric data (subsection 4.1). Second, we also estimate its stellar velocity dispersion 
using the Tully—Fisher relation (subsection 4.2). Then, we estimate the magnification factor 
based on the singular isothermal sphere model for gravitational lensing (subsection 4.3). 

4-1- Photometric Redshift of the Lens Galaxy 

Using optical B , Rc, Ic, NB816, and z' data, we estimated a probable photometric 
redshift of the lens galaxy. Its optical magnitudes obtained in our observations were B = 
25.1 ± 0.06, R c = 24.5 ± 0.04, I c = 24.3 ± 0.06, NB816 = 24.2 ± 0.04, and zJ = 23.9 ± 0.09; 
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the quoted errors are la sky noises. It is noted again that because the lens galaxy is not a 
5-dropout, we had a constraint of z < 2.62. 

A photometric redshift of a galaxy was evaluated from the global shape of the spectral 
energy distribution (SED) of a galaxy, as a redshift with the maximum likelihood, 


L{z,t) = I[exp -- 


2=1 


fi-AFi(z,t ) 


CTi 


( 1 ) 


where /,;, Fi(z,t), and a, are the observed flux, the template flux, and the error of the i-th 
band, respectively, and A is defined as 
EFJi/af 


A = 


E FfM 


( 2 ) 


The SED of a galaxy that we observed was mainly determined by the following four 
factors: (1) the radiation from stars in the galaxy, (2) the extinction by dust in the galaxy, 
itself, (3) the redshift of the galaxy, and (4) the absorption by the intergalactic neutral hydrogen 
between the galaxy and us. We treated the above factors in the following way. 

We used the population synthesis model GISSEL96, which is a revised version of Bruzual 
and Chariot (1993) (see also Leitherer et al. 1996), to calculate an SED from stellar components. 
The SED of the galaxy was determined by its star-formation history. The SEDs of local galaxies 
were well reproduced by models whose star-formation rate declines exponentially (r model); i.e., 
SFR(t) oc exp(— t/r), where t is the age of the galaxy and r is the time scale of star formation. 
Different combinations of t and r generate a similar shape of SED. In this work, we therefore 
used r = 1 Gyr models with the Salpeter’s initial mass function (the power index of x = 1.35 and 
the stellar mass range of 0.01 < m/M© < 125) to derive various SED types. We adopted the solar 
metallicity, Z = 0.02. For these models, we calculated SEDs with ages of t — 0.1, 0.5, 1, 2, 3, 4, 
and 8 Gyr, and called them SED1, SED2, SED3, SED4, SED5, SED6, and SED7, respectively; 
note that the SED templates derived by Coleman et al. (1980), elliptical galaxies (the bulges 
of M 31 and M 81), Sbc, Scd, and Irr correspond to those of SED7 (t — 8 Gyr), SED6 (t = 4 
Gyr), SED5 (t — 3 Gyr), and SED4 (t — 2 Gyr), respectively. We adopted the dust extinction 
curve for starburst galaxies determined by Calzetti et al. (2000) with visual extinctions of Ay 
= 0, 0.2, 0.4, 0.6, 0.8, and 1.0. As for the absorption by intergalactic neutral hydrogen, we 
used the average optical depth derived by Madau et al. (1996). However, Scott, Bechtold, 
and Dobrzycki (2000) showed that the observed continuum depression between 1050 (1 + z) 
A and 1170 (1 + z) A due to the L a clouds (Da) for quasars with z < 3 is lower than that 
expected from extrapolation using the cosmic transmission derived by Madau et al. (1996). 
Further, they also showed that there is a scatter in Da from quasar to quasar. Therefore, we 
also investigated the following two cases: 1) 0.5 x the Madau et al.’s average optical depth, and 
2) 2 x the Madau et al.’s average optical depth. We then estimated the probable photometric 
redshift of the lens galaxy. In this procedure, we adopted an allowed redshift of between z = 0 
and z — 6 with a redshift bin of A z = 0.02. 
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Our results are shown in figure 3. In this figure, we show the distributions of the 
likelihood for the seven SEDs. Our best estimate of the photometric redshift of the lens galaxy 
is %>hot ~ 2.1. Since the likelihood, L(z,T), exceeds 0.1 for z ~ 1.5 — 2.5, it is suggested that 
the lens galaxy is located at 1.5 < z < 2.5. This range is consistent with the information from 
our optical spectroscopy (section 3). SED1 appears to be inappropriate for the SED of the lens 
galaxy. It is found that SED7 also has a maximum likelihood larger than 0.1 around z ~ 3. 
However, since the age of the galaxy at z ~ 3 is considered to be younger than 2 Gyr, we also 
threw out SED7. The luminosity of the galaxy was evaluated to be the normalization factor, 
A, times the luminosity of the model. The range of the evaluated absolute H-band magnitude 
for SED2 through SED6 was between z = 1.5 and 2.5 is —20.4 > Mb > —23.1. 

4-2. Stellar Velocity Dispersion of the Lens Galaxy 

Although it is necessary for the estimate of the magnification factor by gravitational 
lensing to know the stellar velocity dispersion (a v ) of the lens galaxy, there is no direct mea¬ 
surement. Given the redshift of a lens galaxy, we could estimate its absolute magnitude and 
then estimate a probable value of the stellar velocity dispersion. Since it is considered that the 
lens galaxy is located at 1.5 < z < 2.5, we could estimate the stellar velocity dispersion as a 
function of both the redshift and SED based on the Tully—Fisher relation established for a 
sample of galaxies in the local universe. 

Here, we used the Tully—Fisher relation derived by Sakai et al. (2000), 

M b = -8.07(logW 20 - 2.5) -19.88, (3) 

where Mb is the absolute B magnitude and W 2 0 is the full width at 20% of the maximum ve¬ 
locity. Assuming that a v = V mt /\/2 = W 20 / 2 V 2 where V rot is the rotation velocity, we estimated 
a v as a function of z and SED. We evaluated M B using the method given in subsection 4.1. 
The results are shown in figure 4. Since the cases of SED1 and SED7 are rejected in subsection 
4.1, the velocity dispersion is estimated to be in a range between 140 km s -1 and 280 km s _1 
for a redshift range of 1.5 <z< 2.5. 


4-3. Magnification Factor by the Gravitational Lensing 


We are now ready to estimate the magnification factor by gravitational lensing. We 
adopt the singular isothermal sphere (SIS) model for simplicity (e.g., Biuney, Merrifield 1998). 
In this model, the magnification factor for the brighter source can be expressed as 


M, = 


e 


0-0e’ 

where 9 is the angle between the lens galaxy and the source (i.e., SDSSp 
in this case), and the 9 e is the Einstein angle, defined as 

2 A 


9k = 47t 


'LS 


Dos’ 


( 4 ) 

J104433.04—012502.2 

( 5 ) 
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where -Dls is the angular diameter distance between the lens galaxy and the source and -Dos 
is that between the observer and the source. In figure 5, we show the magnification factor as a 
function of z for eight cases of a v = 140, 160, 180, 200, 220, 240, 260, and 280 km s -1 . In this 
figure, we also show the magnification factor for the counter image [M_ = 9/(6e — 0)\, which 
can be seen at 9 = 9^ — (3, where [3 is the angle between the lens galaxy and the intrinsic source. 
If the lens galaxy is located at z ~ 2 and has a v ~ 280 km s -1 , 9e becomes 1"9 and we thus 
obtain M + = infinity. However, if this is the case, since /3 ~ 0 and M_ ~ M + , we would observe 
the counter image with nearly the same brightness as that of the lensed image. Further, if 
M + > 10, M_ is also as large as M + . As we shown in figure 1 that there is no apparent object, 
except the quasar (z' = 19.2 mag) and the lens galaxy (z' = 23.9). Since there is no apparent 
counter image, the magnification factor must be much less than 10. If a v 200 km s -1 , the 
condition of M_ -C M + is achieved at z ~ 2. In this case, we obtain M + ~ 2. 

4-4- Concluding Remarks 

We have found a candidate of a lensing galaxy at 1"9 southwest of the SDSS liigh-z 
quasar, SDSSp J104433.04—012502.2. From our observational data, we find that the lens galaxy 
is located at 1.5 < z < 2.5; the most probable photometric redshift is z p h ot ~ 2.1. Adopting the 
SIS model for gravitational lensing, we estimate the magnification factor to be M + ~ 2 if the 
lens galaxy is located at z ~ 2 and has a stellar velocity dispersion of a v ~ 200 km s -1 . 

WL02 suggested that the SDSS z ~ 6 quasars are systematically brighter than expected. 
Namely, they investigated the probability that a quasar with a rest-frame 1450 A magnitude, 
M 1450 (AB), and a redshift of z was selected into the survey using the survey selection function 
given in Fan et al. (2001). They found that the cumulative probability of M 44 5 o(AB) for the 
z ~ 6 quasars is different from the expected one at the 95% significance level. However, they 
also suggested that this inconsistence is removed if the flux of one or more of the z ~ 6 quasars 
is overestimated by a factor of 2. Therefore, our finding presented in this paper may remove 
the inconsistency raised by WL02. 

It is known that the lensing optical depth increases with increasing the redshift, and 
thus objects at higher redshift are more affected by gravitational lensing from a statistical 
point of view (e.g., Turner 1991; Barkana, Loeb 2000; WL02). As proposed by WL02, it will 
be very important to carefully search for lensing galaxies towards all high-z quasars by deep 
and high-resolution imaging. 

We would like to thank the Subaru Telescope staff for their invaluable help. We also 
thank Dr. Hayashino for his invaluable help. This work was financially supported in part 
by the Ministry of Education, Calture, Sports, Science, and Technology (Nos. 10044052 and 
10304013). 
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Fig. 1. Thumb-nail images of SDSSp J104433.04-012502.2 and its neighbor galaxy (upper panel). The 
angular size of the circle in each panel corresponds to 8". A contour plot of a lens galaxy after the quasar 
had been subtracted is shown in the middle panel. The lower panel shows the spectral energy distribution 
of the neighbor galaxy in the magnitude scale. 
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Fig. 2. Optical spectrum of the neighboring galaxy. The hatched regions indicate the positions of the 
strong night sky lines. The sky spectrum is also shown in the lower part of each panel. 
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Fig. 3. Distributions of the likelihood for the seven SEDs shown as a function of the redshift. The arrow 
in the upper part shows the range of the probable redshift of the galaxy estimated by optical spectroscopy. 
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Fig. 4. Distributions of the velocity dispersion of the lens galaxy shown as a function of the redshift for 
seven SEDs. 


10 
























Fig. 5. Distributions of the magnification factor of the source observed at 9 = l."9 (thick lines) and its 
counter image (thin lines) shown as a function of the redshift for the eight cases for the stellar velocity 
dispersion. We assume the magnification of the former one to be M+, except for the case that 9 is smaller 
than 0 e- 
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Table 1 . List of z ~ 6 SDSS quasars. 


Name 

Redshift 

V 

M b * 

SDSSp J104433.04—012502.2 

5.74* 

19.20 

-28.09 

SDSSp J083643.85+005453.3 

5.82 

18.74 

-28.54 

SDSSp J 130608.26+035626.3 

5.99 

19.47 

-27.73 

SDSSp J103027.10+052455.0 

6.28 

20.05 

-27.23 


*The ^-correction was made by using the composite quasar spectrum in Vanden Berk et al. (2001). 

*The discovery redshift was z = 5.8 (Fan et al. 2000). Since, however, the subsequent optical spectroscopic 
observations suggested a bit lower redshift; z = 5.73 (Djorgovski et al. 2001) and z = 5.745 (Goodrich et al. 
2001), we adopt z = 5.74 in this paper for quasar spectrum in Vanden Berk et al. (2001). 


Table 2. Journal of imaging observations. 


Band 

Obs. date (UT) 

Total integ. time (s) 

mi im (AB)* 

fwhmL (") 

B 

2002 February 17 

1680 

27.1 

1.2 

Rc 

2002 February 15, 16 

4800 

26.8 

1.4 

Ic 

2002 February 15, 16 

3360 

26.2 

1.2 

NB816 

2002 February 15 - 17 

36000 

26.6 

0.9 

z' 

2002 February 15, 16 

5160 

25.4 

1.2 


* The limiting magnitude (3a) with a 2" aperture, 
t The full width at half maximum of stellar objects in the final image. 


Table 3. Photometry of a lens galaxy. 


Object 

B 

Rc 

Ic 

NB816 

z' 

Lens galaxy 

25.1 + 0.06 

24.5 + 0.04 

24.3 + 0.06 

24.2 + 0.04 

23.9 + 0.09 
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